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Recent work on local functional theories of critical inhomo- 
geneous fluids and Ising-like magnets has shown them to be a 
potentially exact, or near exact, description of universal finite- 
size effects associated with the excess free-energy and scaling 
of one-point functions in critical thin films. This approach is 
extended to predict the two-point correlation function G in 
critical thin-films with symmetric surface fields in arbitrary 
dimension d. In d = 2 we show there is exact agreement with 
the predictions of conformal invariance for the complete spec- 
trum of correlation lengths as well as the detailed position 
dependence of the asymptotic decay of G. In d = 3 and d > 4 
we present new numerical predictions for the universal finite- 
size correlation length and scaling functions determining the 
structure of G across the thin-film. Highly accurate analyti- 
cal closed form expressions for these universal properties are 
derived in arbitrary dimension. 

PACS numbers: 64.60. Fr, 05.7Q.Jk, 68.35. Rh, 68. 15. -He 



Finite-size scaling theory predicts that the free-energy, 
order-parameter and correlation length of a critical thin- 
film Ising-like magnet (or of a fluid between parallel- 
plates at the bulk critical point) are characterised by 
universal amplitudes and scaling functions that depend 
only on dimensionality d and the qualitative nature of of 
the boundary conditions . These have been the topic 
of widespread theoretical interest j|] with approaches in- 
cluding sophisticated e-expansion treatments as d^A 
and in particular studies that invoke conformal invariance 
in d = 2 1^ . In a recent article Boran and Upton (BU) 
have revisited the problem of critical thin-films using 
a local-functional model similar to one introduced earlier 
by Fisher and Upton . Specifically BU consider critical 
films of thickness L and infinite area A with both like sign 
(-I--I-) and opposite sign (-1 — ) symmetry-breaking surface 
fields and obtain explicit expressions for both the univer- 
sal Casimir amplitude Aab associated with the finite-size 
free-energy and universal scaling function '^ab{x/L) for 
the equilibrium magnetisation profile m{x), for arbitrary 
dimension d. For the (-1 — h) geometry in particular where 
there is less ambiguity in the choice of functional, the 
predictions of BU are in astonishing agreement with pre- 
vious results obtained from the e-expansion and confor- 
mal invariance. For example in d = 2 the local func- 
tional result for \E'++(a;/L) is exact [|| whilst the numer- 
ical value obtained for is in near perfect agreement 
with the conformal result = — 7r/48 ||]. Whilst the 
local functional theory is phcnomenological in nature, the 



agreement with microscopic theory is compelling and the 
model offers the exciting possibility of predicting (near) 
exact results for universal amplitudes and scaling func- 
tions in d — S. 

In the present work we extend the local func- 
tional theory to allow calculation of the two-point con- 
nected correlation function G(ri,r2) = (TO(ri)m(r2)) — 
(TO(ri))(r7i(r2)) for the {++) thin-film. First we discuss 
the d = 2 Ising strip and compare the theory with ex- 
act results obtained from conformal invariance. This is 
a stringent test of the accuracy of the local functional 
theory and our analysis shows an extraordinary level of 
agreement which we find even more striking than the 
successes obtained by BU for and ^'++(a;/L) men- 
tioned above. In particular we derive the precise spec- 
trum of universal correlation lengths as well as the exact 
and highly non-trivial position dependence of the leading- 
order and next-to-leading order asymptotic decay of G. 
We are also able to obtain the exact Casimir amplitude 
improving on the estimate of BU. Secondly we 
present detailed predictions for the universal finite-size 
correlation length and scaling function "Bdix/L) charac- 
terising the asymptotic decay of G for the d — 3 and d>A 
thin-films. These features of G have not been discussed in 
detail before and we are able to derive accurate analytic 
expressions for them valid in all dimensions. 

To begin we recall some pertinent universal finite-size 
scaling properties of critical thin-films and the exact 
predictions for G in d = 2 arising from conformal in- 
variance. Consider a d dimensional Ising-like thin film 
with surface fields hi, h2 acting on the spins in the sur- 
face planes a; = and x = L. These fields mimic the 
preferential local adsorption of the confining walls on 
a fluid in a parallel-plate geometry. We define the re- 
duced excess free-energy (L) = {F — ALfb)/ AksT by 
subtracting from the total free-energy F the appropri- 
ate contribution from the bulk free-energy density /t,. 
By further subtracting from f^{L) the semi-infinite (re- 
duced) surface tensions of the two-independent surfaces 
we obtain the desired finite-size contribution to the free- 
energy, Af^{l). Finite size scaling theory predicts that 
exactly at the bulk critical point, for asymptotically large 
widths L (and dimension d < 4), this can be written 
Af^{L) ^Aab/L'^~^ + ■ ■ ■ where the Casimir amplitude 
Aab depends only on whether the surface fields have like 
(+-f) or opposite (-1 — ) sign (or are zero). Hereafter we 
restrict attention to the (-I--I-) geometry where the re- 
sults of the local functional theory are most striking. 
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A second universal quantity is the normalised scaling 
function describing the magnetisation profile 

m{x) across the thin-film. For a semi-infinite system the 
asymptotic large distance decay of the critical magnetisa- 
tion profile is m {x)~Cx~^/'' where C is a non-universal 
amplitude and (3, v are standard bulk critical exponents 
for the magnetisation and correlation length respectively. 
In the (-1--!-) critical thin-film finite-size scaling pre- 
dicts that for large x, L with x/L arbitrary the pro- 
file can be written m{x) =m{L/2)'9^+{x/L) where the 
mid-point magnetisation scales as m{L/2) (x L~^^^ and 
^++{x/L) is normalised such that ^'++(1/2) = 1. Fi- 
nally, finite-size scaling theory ||] also predicts that the 
asymptotic decay of correlations along the critical thin- 
film is described by a (true) correlation length with 
a universal critical amplitude ratio. 
It is well known that in d = 2 these predictions are 
beautifully confirmed by conformal invariance and are 
consistent with results obtained from exactly solvable 
models. In particular by using the standard logarith- 
mic function to conformally map the semi-infinite plane 
to the finite-size strip the magnetisation follows [g| as 
m(x) = C(§ sin ^)~^/^ implying ^'++(z) = (sinTra;)"^/*. 
Applying the same mapping to the semi-infinite two- 
point function calculated by Burkhardt and Xue [|l0| we 
obtain 
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coshTry/i — cos7r(a;2 + xi)/L 
coshTry/i — cos7r(a::2 — xi)/L 
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and y denotes the scalar distance between the spins mea- 
sured along the strip. From the asymptotic y/L—^oo ex- 
pansion of (2) it follows that the true correlation length 
in the (H — h) strip is = L/2n identifying the univer- 
sal critical amplitude ratio predicted by finite-size scaling 
theory. More generally, this expansion identifies the spec- 
trum of correlation lengths ^^"^ determining the decay of 
the exponential terms as 



(n) 



L 



{n + 2)7r 
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Thus the explicit conformal invariance prediction for 
the leading-order and next-to-leading-order asymptotic 
y/L—^oo decay of G^^ is 

G++(xi,X2;2/) = — (-1 [sm — sm—j x 



1 + 4cos— -I-cos— -^e ""^'-^ 
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which shows a detailed position dependence across the 
strip. We wish to compare these exact predictions with 
the results of the BU local functional theory. 



The desirable scaling properties of quite general local 
functionals of the magnetisation have been considered 
by Fisher and Upton 0] but in application to the (++) 
thin-film (where there is no zero of the equilibrium pro- 
file) the appropriate model is simplicity itself. BU as- 
sume translational invariance of the trial profiles along 
the thin-film and base their analysis on the excess (re- 
duced) free-energy functional (ignoring boundary terms) 



dx 



K 
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(4) 



where m = d'm/dx, rj and 5 are standard bulk exponents 
and _ftr, u are parameters that must be regarded as in- 
puts into the theory, although their values do not influ- 
ence the determination of the scaling function \l/++(z) 
nor the structural properties of G. The surface fields en- 
ter through boundary terms but these play no role in de- 
termining the universal quantities of interest and may be 
replaced by fixed boundary conditions m(0) = m(L) = oo 
as appropriate to the scaling limit. Specialising to d = 2 
and using the bulk critical exponents /? = 1/8, v = 1, 
77= 1/4 and 5= lb the Euler-Lagrange equation for the 
equilibrium profile is 
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which has a first integral 

Km^^m^ um^^ 



16 



- Ap 



(6) 



where the constant of integration Ap = df^/dL 
—A^+L^'^. From (5) it is straight forward to derive [[ 
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in exact accordance with conformal invariance. From the 
profile it follows that Ap = 7r^if/128i^ implying A^+ — 
— 7r^A7l28. Using scaling arguments BU relate K to 
known bulk critical exponent and amplitude ratios and 
yield a numerical value for extremely close to the 
exact result = — 7r/48 |^ and with an error that 

may even be attributed to uncertainties in the values 
of bulk critical amplitude ratios used to determine K. 
We shall determine K another way and show how the 
exact result for A^j^ may be obtained within the local 
functional theory. 

With these preliminaries over we can now turn to the 
main part of our work and describe the calculation of 
G++ within the local functional theory. As it stands the 
BU functional (4) is not capable of yielding the full cor- 
relation function because it assumes translational invari- 
ance of the magnetisation along the strip. However it is 
reasonable to argue from isotropy that the full functional 
is given by 
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dy dx 
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and recall that we have adopted fixed boundary condi- 
tions as appropriate to the scaling limit. Minimisation of 
(8) leads to the extended Euler-Lagrange equation 



Km ^\7^m 



Km '^(Vm)^ = um^ 



(9) 



which of course reduces to (5) as soon as we impose trans- 
lational invariance along the strip. We mention here that 
the solution to (9) for the magnetisation profile interior 
to a disk of radius R with fixed up-spins at the boundary 
is m(r) = C(f (1 - r2/i?2))-i/8 which is identical to the 
exact conformal invariance result The same is also 
true in all higher dimensions d if we adopt the appropri- 
ate values for the bulk critical exponents in (9). 

Within the framework of density-functional theory the 
two-point correlation function follows from solution of 
the Ornstein-Zernike equation pj| 



5^T[m] 



■G(r2,r3)=5(ri-r2). (10) 



(5m(ri) 5m{Ys) 
For our functional this reduces to the differential equation 
^135 
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G++{xi,x2;y)^5{xi-X2)5{y) (11) 



which we seek to solve using the spectral expansion 



G++(a;i,X2;2/) = V ■(/'«(a;i)C(a;2)e 



-E„y 
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with energy levels En equivalent to the inverse correlation 
lengths l/^("^ discussed earlier. Orthonormality implies 



2K j dx'ipl^{x)m{x) Tpmix) 
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En 



(13) 



which is needed to determine the correct normalisation 
factors of the eigenstates. The spectrum is thus obtained 
through solution of the eigenvalue problem 
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subject to the boundary condition V'n(O) = ipniL) = 0. 
Substituting for m{x) and the value of Ap quoted earlier 
(without having to fix K) this reduces to 



-C(a;)- ^ (cot^) V;(a;) + 
1357r2 



64L2 
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^„(x) = E^^nix). (15) 



which can be solved using standard methods [|^. For 
the eigenvalue spectrum we find 



E„ 



{n + 2)7r 
I 



n = 0, 1, 2, 3, 



(16) 



whilst the eigenvectors can be written as ipnix) 
7V„(sin^)i5/8^„(a;) with 



n n 1 r, nx 
-,2 -f -, -;cos^ — 
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for the even states and 



0„(x) = COS— 
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for the odd states. Here F[a, /3,j; z] denotes the usual 
hypergeometric function and iV„ is an appropriate nor- 
malisation factor which can be found using (13). The 
eigenvalue spectrum (16) is identical to the result ob- 
tained from conformal invariance (2) and exactly iden- 
tifies all the universal critical amplitude ratios for the 
correlation lengths. From the above the asymptotic de- 
cay of the correlation function within the local functional 
theory is given by 



G'^~^{xi,X2;y) = 



C 



sm 



, sm sm I X (19) 
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where C^{K/8uY^'^^ is the non-universal constant for 
the magnetisation while from (13) we find Ni/Nq = 2 in 
precise agreement with the conformal result (3)! More- 
over it is clear that the parameter K simply sets the 
overall scale of the correlation function and that va d = 2 
the appropriate value is if = 37r/8. Utilising this iden- 
tification we are led to the prediction — — 7r/48 
which is identical to the exact conformal invariance re- 
sult We comment here that only small discrepan- 
cies with the conformal result are found at higher or- 
ders which are surely of no practical importance. It is 
also possible to calculate the moments G'J^(a;i, 0:2) = 
j dyy'^'^G^^{xi,X2\y) using the local functional the- 
ory and compare these with the expressions obtained 
from conformal invariance. Consider for example the ze- 
roth moment Gq{xi,X2) from which we can construct a 
universal two-point scaling- function Kq'^{xi/L,X2/L) = 
G'^+{xi,X2)/G^+{L/2,L/2). We omit the details of this 
calculation and simply quote our final result calculated 
within the BU functional [|12| 
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which turns out to be identical to the conformally in- 
variant result obtained by (numerically) integrating (2) 
provided that xi—X2- There are small discrepancies be- 
tween the expressions if xi ^ X2 but even here the local 
functional theory is an excellent qualitative description 
of the structure of correlations across the strip. Never- 
theless this implies that the exact conformal result for 
Gq^ {xi,X2) does not satisfy the correlation function al- 
gebra necessitated by all local functional theories [p^ . 
Given that it is therefore not possible to construct a 
truly exact local functional theory consistent with confor- 
mal invariance predictions for G'++(a:i, X2; y), the level of 
agreement achieved by the theory is all the more striking. 

Next we turn our attention to the d dimensional thin- 
film and focus on the decay of G++(a;i, 0:2; y) for asymp- 
totic parallel displacements |y| 00. Within the local 
functional theory this behaves as 

G-(..„.2,y)-C^s.(^)s.(f)^;;^ (21) 

and recall that we anticipate that is universal. 

Here we have introduced a universal scaling function 
Sj^{x/L) describing the position dependence across the 
thin- film which we normalise (analogous to ^++) such 
that Ed{l/2) = 1. Both Ed{x/L) and ^^^"^ are determined 
by an eigenvalue equation analogous to (15) which de- 
pends only on d, P/iy and 5'++(a;/L). No information 
concerning K or u is required. We omit details of the 
calculation which are obtained assuming the same nu- 
merical values for the critical exponents /3 = 0.328 and 
= 0.632 in d — 3 as used by BU and the standard mean- 
field results (3 = v =1/2 for d>A. For the universal cor- 
relation length critical amplitude we obtain -L/^'-"-' ~7.34 
and L/S}^'^ ~ 8.40 in d = 3 and d> A respectively and 
recall = 27r in d = 2. The numerical results for 

the scaling functions 'E.d{x/L) are quite remarkable and 
are shown in Fig. 1. To extremely high precision (see 
inset) and for all practical purposes these functions are 
completely indistinguishable from the analytic expression 

S,(x/L) = (sin^) (22) 

where d* =min{(i, 4}. Notice this is exact 'vt\ d = 2 and 
embodies the correct short-distance expansion behaviour 
close to each wall Vd. Using 2d(a:/L) we can derive an 
extremely accurate expression for the ground state eigen- 
value and hence obtain for the universal ratio 

eW-v'''* + [^'-^-^ (23) 

where J++(l) = du (1 — u'' )^i/2 is the elliptic function 
introduced by BU. This is exact in d = 2 and within 
0.3% of the numerical values quoted above for (i = 3 and 
d>4. For all intents and purposes this may be regarded 



as the exact analytical local functional prediction for the 
universal critical amplitude ratio. 

In summary we have established that the BU local 
functional model is an excellent description of correla- 
tions in the d = 2 critical thin-film and used the 
theory to make predictions for the universal structure of 
the correlation function in higher dimensional systems. 

EDM and CR acknowledge support from the EPSRC 
and EC (ERBFM-BICT983229) respectively. 
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FIG. 1. Comparison of the universal functions 'Zd{x/L) in 
the d = 2, 3 and 4 (+-f ) thin-films. Inset is shown the discrep- 
ancy A'Ed{x/L) = 'E.d(x/L) — 'Ed(x/L) from the closed form 
approximation (22) which is not visible on the scale of the 
main figure. 
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